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Abstract. We prove that for the two-dimensional steady complete compress- 
ible Euler system, with given uniform upcoming supersonic flows, the following 
three fundamental flow patterns (special solutions) in gas dynamics involving 
transonic shocks are all unique in the class of piecewise smooth functions, 
under appropriate conditions on the downstream subsonic flows: (i) the normal 
transonic shocks in a straight duct with finite or infinite length, after fixing 
a point the shock-front passing through; (ii) the oblique transonic shocks at- 
tached to an infinite wedge; (iii) a fiat Mach configuration containing one 
supersonic shock, two transonic shocks, and a contact discontinuity, after fix- 
ing the point the four discontinuities intersect. These special solutions are 
constructed traditionally under the assumption that they are piecewise con- 
stant, and they have played important roles in the studies of mathematical gas 
dynamics. Our results show that the assumption of piecewise constant can be 
replaced by some more weaker assumptions on the downstream subsonic fiows, 
which are sufficient to uniquely determine these special solutions. 

Mathematically, these are uniqueness results on solutions of free bound- 
ary problems of a quasi-linear system of elliptic-hyperbolic composite-mixed 
type in bounded or unbounded planar domains, without any assumptions on 
smallness. The proof relies on an elliptic system of pressure p and the tangent 
of the flow angle w = v/u obtained by decomposition of the Euler system in 
Lagrangian coordinates, and a newly developed method for the L°° estimate 
that is independent of the free boundaries, by combining the maximum prin- 
ciples of elliptic equations, and careful analysis of shock polar applied on the 
(maybe curved) shock-fronts. 



1. Introduction 

In discussing steady flows which are supersonic in the entrance section of a duct 
and then become subsonic, R. Courant and K. O. Friedrichs wrote in their classical 
monograph Supersonic Flow and Shock Waves [15l p. 372] that 

" We know in general that this assumption (the flow is con- 

tinuous) is not tenable and that we must consider flows involving 
shocks. It is a question of great importance to know under what cir- 
cumstances a steady flow involving shocks is uniquely determined 
by the boundary conditions at the entrance, and when further con- 
ditions at the exit are appropriate." 
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This paper is exactly devoted to showing that, for the two-dimensional steady 
complete compressible Euler system, with given uniform upstream supersonic flows, 
under certain reasonable conditions of the downstream of the subsonic flows, in the 
class of piecewise functions, the uniqueness of the following three important flow 
patterns in gas dynamics: 

(i) The transonic normal shock in a two-dimensional straight duct with finite 
or infinite length, after fixing a point the shock-front passing through (see 
Figure 1); 

(ii) The transonic oblique shocks attached to an infinite wedge against the uni- 
form supersonic flow (see Figure 2); 

(iii) A flat Mach configuration involving three shocks and a contact discontinu- 
ity, where one shock is supersonic, and another two are transonic, with the 
contact line separating two subsonic regions, after fixing a point O where 
the four discontinuities intersect (see Figure 3). 

These flow patterns have been constructed in [15l p. 311, pp. 332-333] by using shock 
polar, under the assumptions that the shock-fronts (contact line) are straight lines, 
and the subsonic flows behind the shock- fronts are also uniform — i.e., the so- 
lution is piecewise constant (see §2 below). The uniqueness results we obtain in 
this work indicate that under certain conditions of the downstream subsonic flows, 
these assumptions may be relaxed to that the flow fields are only piecewise con- 
tinuously differentiable — then for given uniform upcoming supersonic flows, the 
transonic shock-fronts (and the contact line appeared in the Mach configuration) 
must be straight lines and the subsonic flows behind them must be uniform. From 
mathematical point of view, since the steady Euler system for subsonic flows is of 
elliptic-hyperbolic composite-mixed type, these are all results on uniqueness of so- 
lutions of free boundary value problems of a nonlinear system of mixed type, with 
the transonic shock-fronts and contact line being the free boundaries. 




Figure 1. Normal transonic shock in a straight duct. 

These special while important flow patterns have been the focus of research 
for many years, since the pioneering work on transonic shocks in ducts by G.- 
Q. Chen, M. Feldman [3] for the potential flow equation, and S. Chen [10] for 
the Euler system. In a series of papers [SJ HI [3], G.-Q. Chen and M. Feldman 
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uniform supersonic flow 




Figure 2. Oblique transonic shocks attached to a slim infinite sym- 
metric wedge with open angle 26w and zero attacking angle in uniform 
supersonic flows. 



y 




Figure 3. A flat "direct" Mach configuration. Here 5i is a super- 
sonic shock lying in the quadrant II, 52, iSs are transonic shocks lying in 
quadrant IV, I respectively, and C is a contact line. They are all straight 
lines. The flow fleld is piecewise constant. 

established the stability of transonic shocks in finitely long or infinitely long duct 
with square sections or arbitrary sections, under various conditions at the exit 
of the duct. Z. Xin and H. Yin studied similar problems but employed different 
methods [IB]. For the Euler system, H. Yuan showed in [53] various stability and 
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instability results by using Lagrangian coordinates introduced by S. Chen in [11] and 
characteristic decomposition technique [10]. G.-Q. Chen, J. Chen and M. Feldman 
also introduced the method of stream function in Lagrangian coordinates by which 
the Euler system can be written as a single second order equation to study the 
stability, local uniqueness and asymptotic behavior of transonic shocks in infinite 
ducts for the Euler system. Z. Xin, H. Yin and their collaborators also studied the 
problem of stability of transonic shock in nozzles for the steady Euler system, see 

[23 127]. 

For the study of stability of oblique transonic shocks attached to a wedge in 
supersonic flow, there are papers of S. Chen, B. Fang [131 [E], E. H. Kim |20| . 
and H. Yin, C. Zhou [28]. Both the potential flow equation and the complete 
Euler system were used as models and the authors showed stability under suitable 
downstream conditions. 

For the Mach configurations, one has to use the Euler system since there is a 
contact discontinuity. S. Chen firstly studied the stability of a Mach configuration 
under perturbations of upcoming supersonic fiow, for the steady Euler flow and 
pesudo-steady Euler flow [HI [12]. Later on he and B. Fang also showed the stability 
of a wave pattern of regular reflection-refraction of shocks upon an interface [13) . 
which is similar to a Mach configuration. 

Comparing to stability, the progress on the study of uniqueness is rather slow. 
For potential fiows, G.-Q. Chen and H. Yuan [7, proved the transonic shock in a 
two-dimensional or three-dimensional straight duct is unique modulo a translation. 
L. Liu showed uniqueness of subsonic potential fiow in various domains [22], see also 
[23j . By Proposition 3.1 in 9, , in two-dimensional case, many of the uniqueness 
results on subsonic flows also hold for the Euler system, since under appropriate 
boundary conditions, the flow is in itself irrotational, i.e., governed by the poten- 
tial flow equation. However, to our knowledge, there is no any result on global 
uniqueness of transonic shocks for the complete Euler system before. 

In this paper, we are going to establish the global uniqueness, for the complete 
Euler system, of the three flow patterns we previously mentioned. It turns out that 
the method of maximum principles employed in [7'. for potential flow also works 
for Euler flow, but the point is that one should consider now the extreme values of 
pressure p and slope of the flow angle w = v/u, with u, v the velocity component 
of the flow along x,y axis, respectively. In [29], H. Yuan has shown that p and 
w satisfy a first order elliptic system for subsonic fiow. Then as in |2, one needs 
to show if there are nontrivial extremes of p (or w) on the shock-front 5, there 
would be contradictions by a Hopf-type boundary point lemma. It is here we need 
to utilize some nice properties of the Rankine-Hugoniot conditions — particularly, 
a rather simple curve called p — w shock polar (see Figure 4), which shows the 
relation between p and w behind the shock- front. We recommend [15[ 124] for 
detailed computations of these relations. 

This paper is organized as follows. In §2, we formulate the above mentioned 
physical problems (flow patterns) and state rigorously our results, Theorem 12.11 - 
Theorem 12.61 We also introduce briefly the p — w shock polar for the comfort of 
those readers not familiar with this subject. Then we prove in §3 the uniqueness 
of transonic shocks in straight ducts (Theorem l2.1L Theorem [^7^ . and in §4, the 
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uniqueness of oblique transonic shocks attached to an infinite wedge (Theorem [2?3l- 
Theorem l2.5p . Theorem 12.61 on the uniqueness of the Mach configuration studied 
by S. Chen [TT], is estabhshed in §5. 

2. Shock Polar, Special Solutions and Main Results 

2.1. The Euler System and p — w Shock Polar. The complete compressible 
Euler system governing the steady motion of polytropic gas flow in two-dimensional 
space (x, y) G expresses the conservation of mass, momentum and total energy 
(cf. [la pp. 14-23]): 

d4pu) + dy{pv) = 0, (2.1) 
d^{pu'^+p) + dy{puv) = 0, (2.2) 
d,{puv) + dy{pv^ +p)^Q, (2.3) 

\{u' + v^) + ^=b,. (2.4) 

The last equation is also called as Bernoulli law, with 6o the Bernoulli constant 
which is invariant on each flow trajectory, even across a shock-front. Here p, p, c, u, 
V are the pressure, mass density, sonic speed, velocity component along x and y axis 
respectively, of the flow, and 7 > 1 is the adiabatic exponent. For polytropic gas, the 
state function is p = A{s)p'^, with s the entropy, we then have c = ^Jdpjdp — ^p/ p- 
The Mach number of the flow is defined by M = {v? + v'^)/c. When M < 1 (resp. 
M > 1), the flow is called subsonic flow (resp. supersonic flow). It is well known 
that the system ()2.ip - ()2.4p is symmetric hyperbolic for supersonic flow, while of 
elliptic- hyperbolic composite- mixed type for subsonic flow. As in |29j . we introduce 
w — v/uio denote the slope of the flow angle in case u ^ 0, which, together with 
p, plays a key role in our analysis later. 

Let S = {(a;,y) : a; = /(y)} be a shock-front — that is, a curve where the state 
of the gas is discontinuous and physical entropy condition holds (cf . Deflnition 12.11 
in §2.2 below). It is well known that the following Rankine-Hugoniot conditions 
should hold across S [El p. 11]: 

[pu]- f'{y)[pv]=Q, (2.5) 

[pu'+p]-f{y)[puv]^0, (2.6) 

[puv]~ f{y)[pv^+p\=0. (2.7) 

Here [•] denotes the jump of a quantity across S. For example, if the pressure on the 
left hand side of S (upstream) is poj on the right hand side of S (downstream) is p, 
then we set [p] = p—po- Note that bo in (12. 4p is constant along a flow trajectory even 
across the shock- front, which is the Rankine-Hugoniot condition corresponding to 
conservation of energy. 

For a flxed point Q = {f{y),y) on the shock- front, supposing the uniform state 
(po > 0,uo > 0, Wo = = woiPo > 0) of the gas ahead of the shock-front is 
given0, we consider the four algebraic equations p.4p — p.7p . with five independent 
variables p{Q), u{Q), w{Q), p{Q), which are the state of the flow on the right hand 
side of S, and f'{y)- In TP, pp. 306-309], it has been shown that, the flow angle 
9 = arctanw, horizonal velocity u{Q), and density p{Q), as well as the slope of the 

^The flow angle 9 is the angle between the velocity (11,11) and the x axis. Hence w = tan0. 
^Therefore feo in l|2.4|l is ug/2 -|-7Po/((7 - l)po)- 
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Figure 4. A p — w shock polar. C — {0,po) represents the upcoming 
supersonic flow; A — {0,p'^) represents the strongest normal transonic 
shock; w reaches its maximum at B = (w*,p«), and S = {w sonic, Psonic) 
represents a shock that the state behind the shock-front is exactly sonic. 
B' = {—w,,pt) and S' = {—Wsonic, Psonic) are reflections of B,S with 
respect to the p axis. The open arc SAS' represents all possible transonic 
shocks. The open arc CS and CS' represent all possible supersonic 
shocks. The dashed line below C corresponds to those states that do 
not satisfy physical entropy condition. E = {w,pmcak) represents the 
weak shock; F — { W, pstrong ) represents the strong shock. 



shock- front f'{y), can all be expressed by p{Q). Particularly we have the following 
important relation on p and w [151 p. 347]: 



w = ±- 



JL 

po 



7-1 
7+1 



(2.8) 



Here Mq — uq/cq is the Mach number of the flow ahead of 5, and cq ~ vTPo/po- 
Figure 4 shows the curve represented by this expression in w — p plane, which is 
traditionally called a, p~w shock polar. It is easy to see that this curve is symmetric 
with respect to the line w ~ due to "±" in (|2.8p . 

There are several critical points on the p—w shock polar. A is the point where p 
reaches its maximum (denoted by p"*"); B is the point where w reaches its maximum 
w,. The point S = {w sonic, Psonic) has the special property that the state of the 
gas behind the shock- front is exactly sonic, i.e., M ~ 1. For all points on the shock 
polar with p > Psonic (i-e., on the open arc SAS'), the state behind the shock-front 
they represent are all subsonic. B' , S' are the mirror images of B, S with respect 
to the Hne w = 0. 
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For convenience of later reference, we also write out the following relations on 
the state of the flow behind the shock-front [241 P-44, p. 48]: 

u = uo , (2.9) 

Pouo 

(7 + l)p+ (7- l)Po ^2 10) 

sina = — + cota^m (2.11) 

Here a is the angle between the tangent of the shock-front and the x axis. 

2.2. The Physical Problems and Main Results. Now we formulate the physi- 
cal problems we are concerned with, construct the special solutions (flow patterns) 
by using the shock polar, and then state the results on their uniqueness. 

2.2.1. Problem A: Transonic Shocks in a Finite Duct. Let T) — {{x, y) G : —1 < 
X < 1, < y < 1} be a straight duct with finite length. We set = {(— l,y) : 
< y < 1} and Si = {(1,J/) : < y < 1} as the entry and exit of the duct, 
Fq — {{x, 0) : — 1 < a; < 1} and Ti = {{x,l) : —1 < x < 1} as the lower and upper 
wall of the duct, respectively. 

On the walls, it is natural to impose the impenetrability or slip condition v = 0. 
On the entry S-i, we assume the flow is supersonic and moves in V, so the Euler 
system (|2.1|) - (|2.4|) is hyperbolic in the positive x direction. Therefore we impose 
the initial condition p = po > 0,u = uq > 0,v = vq = 0, p = po > 0, with pq, uq, po 
all being constants and satisfying Mq > 1. On the exit, one Dirichlet condition on 
p is given: 

p = Pi, pi is a constant. (2-12) 
It is indicated in |15i p. 373, p. 385] that this restriction on pressure is more physically 



interesting for nozzle flows. Therefore we get a boundary value problem of (j2.ip — 
p.4p in v. We call this Problem A. We are interested in those solutions of problems 
like this with the following structure, i.e., transonic shock solutions: 

Definition 2.1. Suppose 2? is a domain in R^. For a function x = f{y) defined 
on an interval X of R, let 

S - {(/(y),y)e2? : yei}, 

V- = {{x,y)eV : x<fiy)}, 
2?+ = {{x,y)eV : x>fiy)}. 

Then p,u,v, p e L^^{V)r\C^{V^ U 5)nC^(X'+U5) is a transonic shock solution oi 
(|2.iP " (|2.4p if it is supersonic in and subsonic in 2?+, satisfies equations (|2.ip - 
(|2.4p in D~ U and the boundary conditions pointwise in the classical sense, 
the Rankine-Hugoniot jump conditions (|2.5|) - (|2.7p . as well as the physical entropy 
condition across S: 

[p] > 0. 

Remark 2.1. By (|2.7p and entropy condition, there holds f'{y) = puv/ (pw ^+[p]) and 
the denominator never vanishes. Soifp,u,w,p g L^|,(I?)nC^(2?~UiS)nC^(2?+UiS), 
we naturally have f E in Definition 12. II 
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Remark 2.2. For nozzle flows, we know the transonic shock- front will meet the walls 
perpendicularly. So as explained in |291 p. 1348], at the corner points where they 
intersect, the solution will in general be singular. However, for uniform upcoming 
supersonic flow and straight duct, it can be shown that the first order compatibility 
condition holds (cf. Lemma 3.3 in [25l p. 19]), that guarantees the solution u, v,p, p 
to be at the corners, and hence up to boundary in V. For oblique shocks 
attached to a wedge, in our proof we do not require the solution u,v,p,p to be 
up to the vertex, where the shock-front meet the wedge (continuous is enough, 
see Remark 14.11) . Hence to prove uniqueness, our assumptions on regularity in 
Definition 12.11 are reasonable. 

As a simple application of the shock polar, we show existence of special normal 
transonic shock solutions to Problem A. For given uniform upcoming supersonic 
flow (pq > 0,uo > 0,wo — 0,/9o > 0), a p ~ w shock polar is determined (see 
Figure 4), and for w — 0, we have a unique p"*" > pq (i.e., the point A in Figure 4) 
representing a subsonic state, and hence and /' = may be calculated by 

((2^(l2TT0))(l2lT|) with p there replaced by p+. Suppose the shock- front passes the 
point {t, 0) with — 1 < t < 1, then we may take the shock-front a.s x = t, and the 
state behind it is — (j3+, w+ — 0, which is uniform and subsonic. It is 
easy to see that what we get in this way is a transonic shock solution to Problem 
A, with pi = p+ in (|2.12p . We denote this solution by Ub- 

The following theorem is one of the main results of this paper. 

Theorem 2.1. For given supersonic initial data Uq = {po > 0, mq > 0, wq — 0, po > 
0) on S-i, supposing the flow behind the shock-front satisfies u > 0, then Problem 
A has a transonic shock solution {in the sense of Definition \2.1\i if and only if 
pi — p^ in (|2.12[) . // this holds, let U be a transonic shock solution to Problem A 
with shock-front passing through the point {t, 0), then U must be the special solution 
Ub constructed above. 

2.2.2. Problem A' : Transonic Shocks in an Infinite Duct. This is to consider tran- 
sonic shocks in an infinitely long straight duct V = {{x,y) G : — 1 < x < 
00,0 < y < 1}. The boundary conditions are the same as in Problem A, except 
that (|2.12|) is dropped. We call this boundary value problem of Euler system in the 
unbounded domain V as Problem A' . Of course Ub is a solution to this problem. 
Results like Theorem 12.11 also hold: 

Theorem 2.2. For given supersonic initial data Uq = (po > 0, uq > 0,wo — 
0,Po > 0) on S-i, supposing the flow behind the shock-front satisfies u > u^, and 
M < M* with fixed numbers > and M* < 1, then any transonic shock solution 
to Problem A' must be the special solution Ub, provided that the shock-front passes 
the point {t, 0) G V' . 

Like Theorem 1.3 in [7J p. 568], here we do not need any asymptotic condition of 
the flow at inflnity, under the rather strong assumptions that the flow is uniformly 
subsonic behind the shock-front, i.e., M < < 1, and moves uniformly to the 
exit at infinity. 

We also note that Theorem l2.1l and Theorem l2 . 2l implv particularly the instability 
of the special normal transonic shocks under any perturbations of the downstream 
pressure, as have been shown in [71 [29] . 
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2.2.3. Problem B: Transonic Shocks Attached to an Infinite Wedge. Next we eon- 
sider the transonic oblique shocks attached to an infinite wedge against a uniform 
supersonic flow. By symmetry with respect to the x axis (see Figure 2), we may 
consider only the "half-wedge" (a ramp) with an open angle 9w > and zero at- 
tacking angle against the uniform supersonic flow. Let 2? — {{x,y) G : x > 
0,j/ > xtan^v^}, and set Eq = {{x,y) : a; = 0, y > 0} be the "entry" where 
the uniform supersonic flow (po > 0,mo > 0, = 0, po > 0) flows from left 
to right to the ramp W = {{x,y) € : a; > 0,0 < y < xtanOw} with wall 
r — {{x,y) : X > 0,y — xtanOw}- On the wall we have the slip condition 
w = tan 9w- We would like to find transonic shock in the fiow field. We call this 
Problem B. 

Assuming the transonic shock solution is piecewise constant — that is, the shock- 
front is a straight line issuing from the origin O, which is the vertex of the ramp, 
and the subsonic flow behind the shock-front is also uniform, we may obtain special 
transonic shock solution to Problem B by the shock polar. Indeed, for tanOw G 
(0, Wsonic) (see Figure 4), the line w = tanOw meets the p — w shock polar at three 
different points D,E,F. The point D represents a state that does not satisfy the 
physical entropy condition, so we ignore it. The point E corresponds to a weaker 
supersonic shock, that is, the fiow behind the shock-front is also supersonic. This 
supersonic shock has been intensively studied (see, for example, [SI [HI [HI 130] and 
references therein) , and is not of our concern. The point F represents a transonic 
shock which we are interested in, and we may obtain the slope of the shock- front, the 
density and velocity behind the shock- front by ()2.9p - (|2.1ip . We call this solution 
the strong transonic shock Us- It always corresponds to the point {w — tanOw G 
(0, = Pstrong) ou the uppcr part of the shock polar (i.e., on the closed arc AB). 

For w = tan6]Y e {w sonic, w*), both the weak shock and strong shock are transonic 
shocks, and for w = w,, they coincide. We will denote the weak transonic shock by 
Uy,; it corresponds to the point {w = tan 9w € {w sonic, w^),p = Pweak < Pstrong) on 
the lower part of the shock polar (i.e., on the open arc BS). For w > w^,, detached 
bow shock will appear, and for w = w sonic, the weak shock is sonic. We will not 
consider these two cumbersome cases in this paper. 

For the special strong transonic shock, we have the following uniqueness result. 

Theorem 2.3. For given supersonic initial data Uq = {po > 0, mq > 0, wq = 0, po > 
0) on Eo; supposing the flow behind the shock-front satisfies u > and the following 
asymptotic condition: 

w tanOw as {x,y) ^ oo, (2-13) 

here tanOw G (0,w*], then any transonic shock solution to Problem B with the 
restriction at the vertex of the ramp 

W = tan9w, P = Pstrong at O (2-14) 
must be the strong transonic shock solution Us- 

Under a further rather strong asymptotic condition on the shock-front at infinity, 
namely, 

p < Poo as {x,y) ^ oo on S (2-15) 



■^Note that this is the initial data to the Euler system on Sq. 
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with Poo a constant less than p*, we also have uniqueness of the weak transonic 
shock solution. 

Theorem 2.4. For given supersonic initial data Uq — (po > 0, uo > 0,wo = 
0, po > 0) on Y.0, supposing the flow behind the shock-front satisfies u > 0, and the 
asymptotic conditions ()2.13|) . (|2.15p . then any transonic shock solution to Problem 
B with the restriction at the vertex of the ramp 

w = i&ii6w, P = Pweak at O (2-16) 
must he the weak transonic shock solution Uw Here tan^vy € {w sonic tW*)- 

Here we note that by shock polar, one can always determine that at O either 
(|2.14l) or (|2.16p holds. Therefore our requirements in the Theorems are natural. 
We may replace (|2.13p in Theorem 12.31 by 

P ^ Pstrong as {x,y)-^oo (2-17) 

to guarantee uniqueness of Us , and replace (|2.13p in Theorem 12.41 by 

Pweak as {x,y)->-oo, (2-18) 

as well as dropping (j2.15p to show uniqueness of Uw 

Theorem 2.5. For given supersonic initial data Uq = (po > 0, uq > 0,wo = 
0, Po > 0) on So, supposing the flow behind the shock-front satisfies u > 0, and 
the asymptotic condition (j2.17p {resp. (|2.18p ). then any transonic shock solution to 
Problem B with the restriction (j2.14p [resp. (j2.16l) ) at the vertex of the ramp must 
be the strong [resp. weak) transonic shock solution Us [resp. Uw). 

We recall that, in discussing the shocks attached to a wedge, R. Courant and 
K. O. Friedrichs had written [151 pp.317-318]: "Quite aside from the question of 
stability, the problem of determining which of the possible shocks occurs cannot be 
formulated and answered without taking the boundary conditions at infinity into 

account. If the pressure prescribed there is below an appropriate limit, the 

weak shock occurs in the corner. If, however, the pressure at the downstream end is 
sufficiently high, a strong shock may be needed for adjustment. Under appropriate 
circumstances this strong shock may begin just in the corner and thus, of the two 

possibilities mentioned, the one giving a strong shock may actually occur. All 

statements made here are conjectures so far. While there is little doubt that they 
are in general correct, they should be supported, if possible, by detailed theoretical 
investigation." Out results above, may be considered as a part of such investigation. 

2.2.4. Problem C: Uniqueness of a Mach Configuration. Mach configuration is one 
of the fundamental flow patterns in gas dynamics which is closely connected to 
shock reflection-diffraction phenomena. As explained by R. Courant and K. O. 
Friedrichs in [15l pp. 332-333]: "Three shocks separating three zones of different 

continuous states are impossible. Configurations involving three shock fronts 

must therefore involve an additional discontinuity. The simplest assumption, in 
agreement with many observations, is that there occurs in addition to the three 
shock fronts a single contact discontinuity line." Such a wave structure is called a 
Mach configuration (cf. [11] p.2]). Figure 3 shows a possible Mach configuration, 
where Si is an incident shock separating the two supersonic regions flo and f2i; S3 
is called the reflected shock from the contact discontinuity line V; S2 is called a 
Mach front. The flow in il.2 and ^3 are both subsonic. So Si is a supersonic shock 
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and ^2 , S3 are transonic shocks . Following terminologies in [Tl] , we call ilo U 5i U 1 
the supersonic part of the Mach configuration, and S3Uil3U'DUft2US2 its subsonic 
part. 

For the case that all the discontinuities are straight lines and the flow field is 
piece- wise constant, the existence of such a "flat" Mach configuration showed in 
Figure 3 can be proved by using shock polar. In Figure 5, the point C represents 
the upstream supersonic flow in flQ. Ii corresponds to the flow state in fli. By the 
shock polar R roots at /i, we can determine all the possible states in ^,3; and by 
the polar roots at C, we may determine all the possible states in 512- Since w and 
p should be continuous across the contact line T), the point 

l2,3 = {Wm.Pm) (2.19) 

where the two loops intersect represents the w and p in 0,2,^3- Hence by (j2.9p - 
p.lip . the slope of the lines Si, S2, S3, V and u, v, as well as p, in fli, ^2, ^^3, can 
all be uniquely determined. We recommend [TS] pp. 346-350] or [T, pp. 22-24] for a 
detailed discussion. 




Figure 5. Existence of a fiat Mach configuration by using the p — w 
shock polar. 

It arises naturally the question of uniqueness: under what conditions, can we 
show the flat Mach configuration constructed above is the only possible one in the 
class of piecewise functions? We have the follows. 

Theorem 2.6. Suppose the supersonic part, i.e., the uniform upcoming supersonic 
flow in fio and a straight supersonic shock-front Si, are given (see Figure 6), and 
u > 0, as well as the following asymptotic condition in the subsonic region Vt2,^3: 

p^Prn as (x, j/)->oo, (2.20) 

where p„i is determined in (j2.19l) . For definiteness, we also fix the point O € Si 
where it meets the possibly curved discontinuities S2,S3,'D. Then any weak entropy 

solution {p,w,u,p) e L-(R2)nci(r^)nC'HnT)nC'/„,(ri^\{0})nc/„,(ri;\{0}) of 
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the Euler system must be the flat Mach configuration. Here Cl^^{^\{0}) is the set 
of those functions continuous in Q and continuously differentiable in any compact 
subset ofQ\{0}. (|2.20p may be replaced by 

w ^ Pm as U 9 (x, J/) — )■ oo, (2-21) 

and the same conclusion holds. 




For the four Problems stated above, by the theory of classical solutions of hyper- 
bolic systems [H] , it is easy to show that the supersonic flow ahead of the transonic 
shock-fronts must be uniform (cf. jH p. 216] for a sketch). Particularly, we may 
prescribe supersonic initial data on {x = 0} in Problem C (see Figure 6). To prove 
these uniqueness results, one then only needs to show the state behind the transonic 
shock must be uniform, and hence by shock polar, the shock-front must be straight. 
We will show how these hold and prove Theorem 12. 11 Theorem l2.2l in §3, Theorem 
12.31 - Theorem l2.5l in §4, and Theorem l2.6l in §5, by using a rather unified and newly 
developed approach that combining various maximum principles and analysis of 
the shock polar. 

3. Uniqueness of Transonic Shocks in Ducts 

3.1. Reduction of the Euler System. Supposing 

pu>0 and < M < 1, (3.1) 

the Euler system (|2.1I) - (|2.4I) may be written in the Lagrangian coordinates (^, rj) 
as (see [Mj pp.1351-1356]) 

d^P + >^Rdr,p - Pidr,w = 0, (3.2) 
d^w + l32driP + XrOj^w = 0, (3.3) 
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together with the invariance of entropy along flow trajectories 



a? I ^ ) = 0, (3.4) 



and the Bernoulli law 



lu\l + w') + ^ = bo. (3.5) 

Here ^ = x, and measures the mass flux along two flow trajectories. The flow 
trajectory in Lagrangian coordinates is the line rj = constant. The coefhcients in 
([33) and are 

pc^uw 



/3i = > 0, (3.6) 



_ (M2 - l)c^ 

P2 - — ^-n 2V > 

u(u^ — c^) 

These may be obtained from (4.8)(4.19) in [29, pp. 1355-1356], by taking ip = ip' ^ 
there. Note our assumption < Af < 1 in p.ip implies m 7^ c, which is required 
in (4.21) in [29, p.l356]. By (3.11) or (3.12) in [29, p.l352] and our assumptions 
in Definition 12. 1[ the transformation {x,y) 1— >■ (^,77) is and invertible, provided 
that pu > 0. It is easy to check that for subsonic flows, (|3.2p and (|3.3p constitute an 
elliptic system, which is the elliptic part of the Euler system. Equations (|3.4I) (13. 5p 
are the hyperbolic part of the Euler system. 

The equations (13. 2p p.3|) may also be expressed as 

Dw = WDp, (3.7) 

with Dw — {d^w, drjw)'^ the gradient of a function w, and W the matrix 



Note that detW = P2/P1 = (1 - Af^)/(p^u'*) > 0, so is invertible. By acting 
(9,,, —d{) from left to p.7p . we get a second order elliptic equation of p in divergence 
form: 

+ Y^d.p) + d, + (/32 + ^)a,;,) = 0. (3.9) 

By our assumptions in Definition l2.1[ the coefficients^ are all functions^, and the 
equation is uniformly elliptic in any fixed compact domain. Therefore the strong 
maximum principle for weak solutions of elliptic equations in divergence form is 
available, see Theorem 8.19 in [19l p. 189]. We also need the following boundary 
point lemma due to R. Finn and D. Gilbarg (see Lemma 7 in [HI p.31]). 

Lemma 3.1. Let the coefficients of the system 

— Vy — aUx + h'uy, Vx — b" Ux + cUy, (3.10) 
A = 4ac-(6' + 6")^><5>0 (3.11) 



^They are \r/Pi, P2 + A|j//3i. 

^Hence bounded in any fixed compact domain. 

^Although the formulation here is a little different from that in I18| , this lemma is obviously a 
direct consequence of that. 



14 



BEIXIANG FANG, LI LIU, AND HAIRONG YUAN 



be Holder continuous in the closure of a simply connected domain fl C with 
C^'" boundary (0 < a < 1), and u,v E C^{il) is a solution to p.lOp . Suppose there 
is a Q (z dQ satisfying u{x) > u(Q) {u{x) < u(Q)) for all a; G 57, then there holds 
rill rill 

|^(Q)>0, (|-(Q)<0,) (3.12) 

where n is the inward drawn normal of the curve dCl at Q. 

It is clear that the inner normal n in (I3.12p may be replaced by any vector I at 
Q with lun > 0, since along the tangent r of dfl we have ^{Q) = 0. In this paper 
we always use "•" to denote the scalar product of vectors. 

3.2. Finitely Long Duct: Proof of Theorem 12.11 Under the assumptions in 
Definition 12.11 and Theorem 12. 11 we show (|3.ip holds. We directly have u > and 
M <1. By (I2H), ([2T01) and ([331), the entropy A{s) = p/p^ has the estimate 

< < A{s) < s* in V+ (3.13) 

with two numbers s*, s* depending only on the upcoming supersonic flow pq^uq, po 
and the adiabatic exponent 7 > 1. Hence from u > and M < 1 we get < 



u 



2 < = jA{s)p''~^, therefore p > in 2?+. One then checks that for ([3J1) . 
A = 4/32//?i = 4(1 — M'^)/{p'^u*) is bounded away from zero in a compact subset 
of So we may employ the results in §3.1 for the following analysis. 

The duct, in the Lagrangian coordinates, is now the rectangle {{£,,ri) : —1 < ^ < 
1, < 77 < 770}, with r]o = pqUq > (see Figure 7). At the exit Ei — {(1, 77) : < 
V ^ Vo} we still have the Dirichlet condition for pressure: 

P^Pi, (3.14) 

with pi a constant (cf. (12.121) ). On the upper and lower walls Fi — {(■^,770) '■ ~1 < 
f < l},ro = {(C, 0) : -1 < ^ < 1}, there should hold w = 0, hence d^w = 0, and 
by p.7p , this is a Neumann condition of p: 

driP = 0, (3.15) 

since now we have Xr — Q and ^1^2 > 0. 




Figure 7. The duct and transonic shock in (^,77) coordinates. 

Without loss of generality, we suppose the shock-front S, with the equation 
^ = tpiv) £ C'^[0i'7o] in Lagrangian coordinates, passing through the origin O 
(i.e., 7/^(0) = or /(O) — 0). We do not need write out the Rankine-Hugoniot 
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conditions in Lagrangian coordinates here because the relations (|2.8|) — (|2.10p do 
not depend on coordinate system (these Rankine-Hugoniot conditions are available 
in p. 1358]). Since w = on the walls, the shock- front always meets the walls 
perpendicularly (i.e., ip'i^) = = ip' {tjq)), even in the (^,77) coordinates Q. 

Let n be the domain bounded by the shock-front, walls and exit, where the flow 
is subsonic. Now if p is a constant in fl, that is, p = pi according to p.l4|) . then 
by (|3.7p . is a constant, and by boundary conditions on walls, we have w = 0. 
So from the shock polar, we can infer that pi — p^ . If pi ^ then we get a 
contradiction and either the solution is not constant or C^[Vl) solution does not 
exist. We further may solve on S that p — p'^,u — — p+ and /' = by 

(|2^ — ([2lT|l . and by the hyperbolic part (|01) - ([331) . obtain that p = p+,u = u+ 
in f2. Hence if p is a constant in Q, then pi — p"*" and the solution is exactly the 
special solution Ub- The uniqueness is therefore proved. 

Now we show that if p is not a constant in fl, then there will occur contradictions 
by using maximum principles. In the following, we always assume that p is not a 
constant in Q. 

Step 1. First, by (|3.9p and the strong maximum principle (Theorem 8.19 in [19} 
p. 189]), Pmax, Pmim the maximum and minimum of p in Q, respectively, can only 
be attained on the boundary d^l = Si U (Fq n fi) U (Fi n il) U 5. Since pi is a 
constant, Pmax and Pmm can not be attained simultaneously on Si. By Lemma l3.1l 
and p.lSp . none of them can be attained on the walls. So at least one of them 
should be attained on S. 

Step 2. Now we suppose that Pmax is attained on S. By shock polar, we can 
see the maximum is p^ (corresponding to w = 0, that is the point A in Figure 4) 
and is certainly attained at O = (OjO) and H — {'ilj{rio),rio). Let r — (i/''(?7),l) 
and n — (1, —Tp'^rj)) be respectively the tangential and inner normal vector of S at 
{'ijj{ri),ri), where p attains its maximum. By p.7p . we have 

^ = tDw ^ {tW)Dp. (3.16) 

OT 

But 

[tW) • n = tWt7 = ^((AiiV' - 1)' + /3i/32^") > 0, (3.17) 
Pi 

by Lemma |3. II we get 

{tW)Dp < 0. (3.18) 

Note that at the corner points O, if, we cannot employ Lemma 13.11 directly. 
However, since both the corner angles are 7r/2, we may use a reflection argument to 
overcome this difficulty. Taking O as an example. We first refiect 5, with respect 
to 77 = and obtain its image 5', fi'. Since the corner angle is 7r/2, the curve SVJS' 
is still . Let Br{0) be an open disk centered at O with a very small radius r, 
we set O = Br{0) n (fi U fi' U Fq). Then we use even reflection to extend p (both 
the solution p and those appeared in the coefficients through (|3.7I) ) to O, and by 
Neumann condition (|3.15p . p £ (7^(0). Similarly we use even reflection to extend 
/5, M, odd reflection to extend v,w to O. Hence in O equation p.9p still holds and 



''See (6.5) in [29l p.l358], that is, tp' {r}) = [uw]/\p]. 

^Here the walls are the line segments n f2 (i = 0, 1) excluding their extreme points, hence 
they are open. Note that, however, by our definition, Ei and S are closed. 
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the coefficients are at least Lipsclitiz continuous. We then use Lemma [XT] at O E O 
to obtain (|3?T8l) . 

Hence whenever p reaches its maximum, we have ^ < by (IXTC)) and (IXTa . 
Thus from u;(0) — 0, we may infer that w is always negative on iS\{0}. This is 
impossible since = 0. 

Step 3. Next we consider the case that only the minimum of p is attained on 
S. If Pmin = P'^ , then by shock polar, we see that on S there should always hold 
p = p^ and w = 0, hence ip' ~ 0. We then get a contradiction from p.l6|) and 
p.l7p : the left hand side of p.l6|) is zero, while the right hand side should be 
positive due to Lemma |3. II 

Step 4. We then consider the case that Pmin < ■ First we show Pmin 7^ P*- 
Indeed, if Pmin = P*i by shock polar there are only two possibilities: the points 
corresponding to B or B' in Figure 4, where p = p^,, and w attains its maximum 
w,j(at B), or minimum —w^ (at B'). By (13.171) . we still see the left hand side of 
p. 161) is zero, while the right hand side is positive, a contradiction! 

Step 5. Then we conclude that p„iin G (PsonicP*) U (p*,p'^). We remark that 
Pmin should bc larger than p sonic by Definition 12. 1[ since for those points on shock 
polar below S and S", the flow is supersonic, and sonic at S", S". 

Considering p as a function of w forp S {p sonic, P*) andp £ (p*,p+) respectively, 
we see that p' {w) ^ 0. Now suppose p attains its minimum Pmin at a point P on 
S, then at P we have 

= tDp = t{p'{w)Dw) = p'{w)tDw = p'{w)tWDp (3.19) 

due to (|3.7[) . Therefore rWDp = 0, a contradiction to Lemma [Ol bv (|3.17p . which 
yields rWDp > 0. This finishes the proof of Theorem 12. II 

Remark 3.1. We may replace ()2.12|) or p.l4p . the Dirichlet condition of p, by w = 
on the exit Si. The proof is similar to those presented below in §4.1, and we omit 
the details. 

3.3. Infinitely Long Duct: Proof of Theorem 12.21 For infinitely long duct, 
we set n' as the subsonic region behind the shock-front S: 

^' = {(e, v)-^> ^Kv), 0<V< Vo}- (3.20) 
We first show that under the assumptions in Theorem l2.21 there holds 

\\p,p,u,v\\c^nr)<Co, (3.21) 

with Co depending only on po,uo,pQ and 7. Indeed, by Bernoulli law, the bounds 
ofu,v and c is max{V2, ^7 ^ But p= (0^/(7^(5))) by (|3?T3l) . p is 

bounded. Hence due to p — c^p/j, p is also bounded. 

We claim the coefficients in equation p.9p is bounded, with a bounds determined 
by po,po,uo and 7,m*,M*. To prove this, from p.l3|) . we have a lower bound on 
P- 
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Then by (|3.6|) . it is easy to check that 



1 - M2 - 1 _ M 



■ 2 ■ 



,2„3 / ^ -0 
1 



piK < i3i < 5, 



</?2 < 



Co - u,(l-Af*2)' 

and the claim is proved. 

Finally, since M < M* < 1, it is straightforward to check that the equation p.9p 
is uniformly elliptic in ft' , and p. lip holds in fi'. 

We now state an extreme principle of bounded solutions to general linear elliptic 
equations in the unbounded domain fl' with Neumann conditions on the lateral 
boundaries. 

Lemma 3.2. Let Lu = d^{a^^d^u + a^^d,ju) + djj{a^^d^u + a^^djju) — be a 
linear uniformly elliptic equation with bounded coef ficients a'^ G C/„^(r2')-0 

Suppose 

u S dp.') n Cl^JVl') is a hounded weak solution to this equation that satisfies 
the homogeneous Neumann conditions on the lateral boundaries Fq = {(^, 77) : ^ > 
■0(0), ?7 = 0} and T'l = {(C,?]) : i > il>{rjQ),rj = rja}. Then either Umf — infj^w or 



u 



sup 



supjp-it is attained on S. 



By this lemma, we see at least one of the extremes, Pinf and Psup, which are 
respectively the infimum and supremum of p in fi', is attained on S. The analysis 
in §3.2 then shows the uniqueness of the transonic shock solution, hence Theorem 
HJis proved. B 

To prove Lemma 13.21 we need an up to boundary Harnack inequality for elliptic 
equations with Neumann boundary conditions: 

Lemma 3.3. Let il C R" be a domain with part of its boundary T C lying in 
a plane, and X]"j=i di{a''Wju) = a linear equation with coefficients a*^ satisfying 

n 

^ a^'Uj > A|CP for any (^i, •••,$„) € R", (3.22) 



\a'^ < A in n (3.23) 

j,i=i 

for two fixed positive constants A, A. Suppose u G is a continuous nonneg- 

ative weak solution satisfying the homogeneous Neumann condition on T , that is, 
for any € C^(R") whose support suppi/? is disjoint with dfl\T, there holds 



L 



J2 a'^djud^tpdx = 0. (3.24) 
suppi^n(nuT) ; , , 



^Cf^^{Q') is the set of those functions which arc in any compact subset of fl' . 

1°A W^^^'^in') function u is a weak solution to a mixed boundary value problem if it satisfies 
the Dirichlet condition on S in the sense of trace, the Neumann condition on ro,ri as stated in 
Lemma 13.31 below. 

^^It is worth to note that we need not assume u,v,p,p G C^{Q') in the proof of Theorem 12.21 
The assumption u,v,p,p G ^lotfP''^ °° regularity is enough. 



18 



BEIXIANG FANG, LI LIU, AND HAIRONG YUAN 



Then for any compact subset K of flU T, there holds 

max It < CminM, (3.25) 

K K 

with the positive constant C relying only on K,fl,n,A/X. 

Proof. Let Br{Q) be the open ball centered at Q € R" with radius r. For K ~ 
Br{Q) CC 51, p.25p is the classical interior Harnack inequality (Theorem 8.20 in 

m p-199]). 

We then consider the case that Q E T. Without loss of generality, we may 
suppose T lying in the hyperplane Rq :— {x — {xi, ■ ■ ■ ,Xn) G R" : Xn = 0}. By 
even reflections of a*-'(l < i, j < n — 1 or i ^ j — n) and u with respect to Rq, 
and odd reflections of a*^(i = n or j — n but i ^ j) with respect to Rq, they are 
all extended from B^^{Q) := B2r{Q) n {a; = (xi, • • • , x„) G R" : x„ > 0} C 11 
to the whole ball B2r{Q) for a small r. We also get an elliptic equation satisfying 
p.22p p.23|) in i?2r(Q)- It is straightforward to check that u is still a nonnegative 
weak solution (in the sense of (8.2) in [19l p. 177]) to this equation. Hence by the 
interior Harnack inequality we also have max-g-^^u < C min-g-jgju, and p.25p 

follows with K = Bt{Q). 

For K any compact subset of il U T, we may cover it by finite balls (or half balls 
with center on T) , and p.25p holds by an argument as that in the proof of Theorem 
2.5 in 19, p.l6]. □ 



Proof of Lemma \S. 21 If u is a constant, then obviously Lemma [3.21 holds. In the 
following we always assume that u is not a constant. We prove by contradiction, 
that is, we assume that 

Uinf < minw, Ug^p > maxw. (3.26) 

The argument below is similar to that of proof of Theorem 1.2 in [71 p. 570]. 

Set n'^ = 11' n < L} and Sl ^ n' D ^ L}. By (jX^ . and applying 
Lemma |3. 11 as well as the strong maximum principle in 11'^, we infer that there is 
an increasing sequence {Lk} tends to infinity, and there are points Pk,Qk G ^l^j 
such that 

■u(Pfe) = maxM Usuv, (3.27) 

SI' 

u{Qk) = minu\Mj„/ (3.28) 
n' 

as fc — >■ 00. 

We may choose Li rather large. Then define Uk — {{(,,11) E ft' : Lk — 2 < ^ < 
ifc + 2,0 < ry < r/o}, and Vk = € 11' : Lfc - 1 < C < ifc + 1, < 77 < 770}- By 

translation along ^: 

e - e' + ife, v = v', 

we may translate Uk to Un = {{CiV') : ~2 < ^' < 2, < 77' < 770}, and Vk to 
Vo — {{£.', 1]') ■ < < ^,0 < v' ^ Vo}, and the equation iu = in Uk are 
transformed to L'i.Uk = in [/qi with Ufc(^', 7/) = u(^' + Lk, rj'). 

Now let Vk = Usup — Uk, which is positive and also a weak solution of L'f.v — 
in Uq. Applying Lemma 13.31 to Vk, by the strong maximum principle, we have 

Usup " u{Qk) < maxwfe < Cminvk < C{usup ~ '"(^'a;))- (3.29) 



UNIQUENESS OF TRANSONIC SHOCKS 



19 



It is crucial to note that by our assumptions in Lemma I3.2| the constant C here 
actually does not depend on k. Then taking fc — > oo, we have 

< Usup - Uinf < C X = 0, (3.30) 

a contradiction desired. This finishes the proof of Lemma [3^ 

4. Uniqueness of Transonic Shocks Attached to a Wedge 

4.1. Uniqueness of Strong Transonic Shock: Proof of Theorem [273l Now 

we consider the oblique transonic shocks attached to an infinite wedge against 
uniform supersonic flow (see Figure 2). 




Figure 8. The oblique transonic shock attached to a wedge in La- 
grangian coordinates. 

In Lagrangian coordinates, the wall of the wedge is F = {(^, r/) : ^ > 0,1] — 0}, 
see Figure 8. On F the slip condition is still 

w ~ tan6'w- (4.1) 

As in §3.2, let the equation of the shock-front 5 be ^ = ■ip{i]) with ?7 > and 
-0(0) = 0, and = {(^, ry) : 77 > 0, ^{r]) < £, < 00}. Condition (|2.14|) is now 

W^tanOw, P^Pstrong (4.2) 

which holds at the origin O, and (|2.13p is replaced by 

w^tanOw as (^,?7)->oo. (4.3) 
By the latter requirement and shock polar, we infer that on S, there holds 

P^Pstrong Or p Pnjeak aS {^,7j)->-CO. (4.4) 

We derive a second order elliptic equation of w. Indeed, by p.7p . we have 

Dp = W^^Dw, W-^ = -^W, (4.5) 

P2 

then w satisfies 
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and p.l7p is now 




(4.7) 



One may check that Lemma |3 . 1 1 holds for (j4.5p . and the strong maximum principle 
holds for (gH). 

If It; is a constant, i.e., w = tanOw in 57, then by (|4.2|) and shock polar, we 
infer that there should hold p = Pstrong on S and S is the special strong transonic 
shock- front by (|2.1ip . From (|2.9p and (|2.10l) . u and p can also be calculated on S. 
Then using the hyperbolic part (|3.4l)(|3.5p . we show the solution must be Us, and 
Theorem 12.31 is proved. 

In the following we show there are contradictions if w is not a constant. Suppose 
now that w is not identical to ianOw in Then by Lemma 13.11 and the strong 
maximum principle, either the supremum of w in f2, Wsup, or the infimum of w in 
fl, Winf, must be attained on S. (By (|4.ip and (14.31) . Wsup is indeed the maximum 
of w, while Winf is the minimum oi w.) 

We consider firstly the case that tan6'w G (0, w*). 

Step 1 . Suppose Wsup is attained at a point P on 5. If tan 6w < Wsup < w* , then 
P O. Note that by shock polar, we may regard w as a function of p, w = w{p), 
which shares w'{p) ^ for w ^ ±w*. Now considering the directional derivative of 
w along S at P, ^(P), which vanishes, we have 



Hence tW~^Dw(P) = 0. However, by (jiJ]) and LemmalSJl one gets tW^^Dw{P) > 
0, a contradiction! 

Step 2. Secondly consider the case Wsup = w*, whence p{P) = p*, P ^ O. Then 
by (j4.5p and Lemma IXTl we have 



This means p is strictly increasing once it attains p^. However, let P be the first 
point on S where w attains its maximum as rj varies from to oo. Note that 
Pstrong > P*, SO p dccrcascs near the left side of P, we should have TDp{P) < 0, a 
contradiction to the assumption of p € (5) . 

Step 3. We then conclude that Wsup = tan^iv, which is attained at O. So 
Winf < tan^vi^, and definitely Winf is attained on 5\{0}. If Winf 7^ — and 
Winf 7^ 0, then as in Step 1 above we can show a contradiction since w'{p) 7^ 0. 

Step 4. Now for Winf = — w*, as that of (|4.9p . we may get rDp — tW^^Dw < 0. 
Therefore, whenever w attains its minimum — w*, p is strictly decreasing (hence less 
than p*) and w should run into (— w*, —Wsonic) (the arc B'S' of shock polar), this 
contradicts to our asymptotic condition (j4.3p : w — >■ tan^^K > 0. 

Step 5. For the case Winf — 0, where p attained its maximum on 5\{0}, we 
get = rDp = tW~^Dw < 0, a contradiction. 

Secondly, for the case that tanOw = w*, we have w{0) — Wsup — by (|4.2p 
and (|4.3p . This fact, together with (|4.ip . (|4.3p and strong maximum principle, 
implies Winf < w* should be attained on 5\{0}. Wc then repeat step 3 - step 5 
above to reach contradictions. 

This finishes proof of Theorem 12.31 



tDw{P) 



w'{p)tDp{P) = w'{p)tW-'^Dw{P). 



(4.8) 



rDp = tW-^Dw > 0. 



(4.9) 
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4.2. Uniqueness of Weak Transonic Shock: Proof of Theorem 12.41 Next 
we consider the uniqueness of weak transonic shock. Now at the vertex of the ramp 
O, we have 

W = ta.new, P = Pweak (4.10) 

by (I2.16p . The condition at infinity is (|4.3p plus 

p < Poo as 77 — > cx) on S (4-11) 

with Poo a constant less than . 

As in §4.1, to prove uniqueness, we only need to show that there are contradic- 
tions if w is not a constant, and then either its supremuni Wsup or infimum Winf is 
attained on S. 

Step 1. Now consider the CcLSG Slip IS attained on S. If tawOw < Wsup < w*i 
the proof is the same as in Step 1 in §4.1. For Wsup = w*, by (|4.9p . once w reaches 
its maximum, p should increase (larger than p^). This would be a contradiction to 
our assumption (14.111) and ()4.4p . 

Step 2. For the case that Wmf < ta,n0w is attained on S, if Winf lies on the 
arc BAS' , by continuity, w should first attain its maximum (note now w runs 
counterclockwise on the shock polar from {iSinOw :Pweak) to B' on the arc BAS'). 
There is a contradiction as we showed above in Step 1; if Winf lies on the arc SB^ 
just repeat Step 1 in §4.1 for a contradiction. 

Theorem 12.41 is also proved. 

4.3. Proof of Theorem l2.5l Finally, we show Theorem [23] holds. The arguments 
are very similar to that in §3.2, so we just outline the main ideas. 

By (|3.7p and (|4.ip . one gets an oblique derivative condition of p on W: 

^Rd^P + + \l)d,,p = 0. (4.12) 

The asymptotic conditions p.l7p (|2.18p are respectively 

P~>Pstrong aS (^, 77) CXi, (4-13) 

P^Pweak as (^,77)^00. (4.14) 

To show uniqueness, as before, it is sufhcicnt to show contradictions if p is not 
constant. In the following we always assume p is not identical to Pstrong (resp. 
Pweak)- Then by ()4.12p - ()4.14p and strong maximum principles applied to equation 
p.9p . either Psup — supj^p or Psup = inf^P is attained on S. 

4.3.1. Uniqueness of Strong Transonic Shock. Step 1. Consider the case that 
Psup > Pstrong is attained at a point {Tp{r]^,), 77* ) G 5. If Psup = P^, then as of Step 2 
in §3.2, w should be negative on 5 for 77 > 77*. Note that (|4.13p implies particularly 
w — >■ ±tan6'H'- So we get w — >■ — tan6'iy < 0. Hence V''('7) — uw/[p] has a 
negative upper bound for large rj, which means S will eventually intersect with the 
boundary {(^, 77) : ^ = 0, 77 > 0} where hyperbolic initial data is prescribed, that is 
not allowed. 

Step 2. If p+ > psup > Ps trong, then Psup Hgs On the arc AB oi p — w shock 
polar, and by (j3.19l) there is a contradiction. So Psup = Pstrong- Since p is not 
constant, we infer pinf < Pstrong and is attained on iS\{0}. Just repeat Step 



See footnote 5. Note that [p] > and u > are both bounded away from zero on the 
shock- front which can be checked by using shock polar. 
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4-Step 5 in §3.2 for contradictions. This finishes the proof of Theorem [23] for the 
case of strong transonic shocks. Note that the arguments here also work for the 
case tan 9w = w,, . 

4.3.2. Uniqueness of Weak Transonic Shock. Consider as weh the case that psup > 
Pweak is attained at a point {ip{v*)j V*) ^ ^- If Psup = P'^ , then the argument is the 
same as in Step 1 in §4.3.1. li Psup = P*, refer Step 4 in §3.2 for a contradiction. 
If Psup S {Pweak,P*) U (p*,^"''), refer Step 5 there. Hence we conclude that Psup — 
Pweak, and pinf < Pweak is attained on 5\{0}. Therefore p™/ G iPsonie,P*) and 
refer also Step 5 in §3.2 for a contradiction. This finishes proof of Theorem [23] for 
the case of weak transonic shock. 

Remark 4.1. In the above proofs we do not need to take care of the corner points 
O, so we may only require that p, p,w,u G Cl^^{n\{0}) n C;^^(17)F^ 

5. Uniqueness of a Mach Configuration 

Since the contact line V is also a trajectory of the flow, by virtue of Lagrangian 
coordinates, it may be transformed to the ^ axis, see Figure 9. In the following 
we use superscript -^^^ such as w^'^p^*^, p^'^ to denote quantities in fli, with 
i = 2, 3. I? is a contact line means 

^(2)^^(3)^ p(2)=p(3) onr? = 0. (5.1) 

That is, w and p are continuous in ^2 U fis, although others, such as, p, u, v, s, are 
not. The supersonic initial data is prescribed on the rj axis. 




Figure 9. The possible curved Mach configuration in Lagrangian co- 
ordinates 77). 

We note that by the analysis of shock polar in §2.2.4, we can uniquely solve 
p(2)(0),u;(2)(C)),w(2)(C)),p(2)(0), (O) , ^(3) (O) , ^(3) (o), p(3) (Q), as well as the 



'Cioci^) is the set of those functions that are continuous in any compact subset of f2. 
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slope of the shock- fronts 82,33 at O. Particularly, we have (cf. p.l9|) ) 

pi^){0)=/^HO)=p„., w<^'\0)^w^^\0)=w^. (5.2) 
In the subsonic region — 2,3, the Euler system reads 

a^p^') + Alj^9r,p('' - /3f^9^w(') = 0, (5.3) 

d^w'-''^ + f3^"^dnP^'^ + A^^a^w(*) = 0, (5.4) 

9,((^)=0, (5.5) 

i(«W)2(l + (u,W)2)+ (£^ = 5„. (5.6) 



Here the coefficients are 



pW(cW)2u(i)u,(i) 
(uW)2- (cW)2 ' 

_(pri)2(e(^))2(^W)3 
(wW)2 - (cW)2 

((MW)2 - l)(cW)^ 
uW((uW)2- (cW)2) 



(5.7) 



/3« = 



Now we begin to prove Theorem l2.6l It is divided into several steps. We consider 
only the case of the condition ()2.20p . The case of p.2ip may be proved similarly 
by adopting arguments in §4.1 and then combining the following reasoning, so we 
omit it. 

Step 1 . If both p^^^ and p'^^-' are constants in , ^^3 respectively, then we see the 
subsonic flow in D,2, are uniform, and by shock polar, the transonic shock- fronts 
are straight lines, and then according to analysis in §2.2.4, the Mach configuration 
must be the flat one constructed there. Theorem 12.61 is then proved. 

Step 2. Now we show there will be contradictions if either p^^^ or p'^^^ is not 
constant. There are two cases: I. One is constant while the other is not; II. Both 
are not constant. 

Step 3. Case I. Without loss of generality, we assume p*^^-' is a constant, but 
p'-^^ is not. Then p^^^ = pm in Q2: hence p'^^^ = pm on = 0. 

Let piiip and p^^j: be respectively the supremum and infimum of p^"^-* in fl^. If 
Psup > Pm or p^^^'j: < Pm, by the strong maximum principle, they can only be 

attained on the shock-front 53\{0}. As shown in §4.3, there will be contradictions. 

(2) (2) 

Step 4. Case II. Denote psup and p^^j: be respectively the supremum and 
infimum of p^^^ in Q2 ■ 

(21 

Step 4.1. We firstly show psup > Pm implies contradictions. Suppose now 

(2) . . . (2) 

Psup > Pm- Then by the strong maximum principle, the extreme psup can only be 

attained on the shock-front S2\{0} or on the contact line 77 = 0. As shown in §4.3, 

there will be contradictions if this happens on So it is only possible to take 

place on 77 = 0. By (|5.ip . we infer piup > Pm- As before, pi^p also can only be 

attained on 77 = 0. Hence we get 



p(3) ^ (2) 



(5.8) 
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and they are attained at the same point Q [Q ^ O) on the contact hne. Therefore, 
by Lemma l3.H we have 

V'HQ)<0, (5.9) 

5„p('HQ)>0. (5.10) 

Now let us look at the equations (|5.3p(l5.4p . We may solve (cf. (|4.5p ) 

= ~^{d^w^''> + A(^^9,ii;(')), (5.11) 

= -^(a(^)5,«;« + (/?«/?« + {\ff)d,w(\ (5.12) 

Then by (gSlldOl) and ((O)) (|5TU)) . we may get, at Q, 

(aeu;(3),a^u;(3)), (i,Ag^) > 0, (5.13) 

(a5«;(2),a^«;(2)),(l,A^'^) <0. (5.14) 

On the other hand, note that we have d^p''^\Q) = d^p''^\Q) = by (EH), therefore 
at Q, 

(9^u;«,a,«;«) . (A«,/3«/3« + (A«)2) = 0. (5.15) 

We then get 

d^w^'^Q) = (o»4W«,a^u;(*)).(l,0) 

Pi "2 

Pi P2 



<0 i = 2. 



2 



(5.16) 

However, by (15. ip . we should have d^w'^^^Q) = d^w'^^^Q), a contradiction! 

Step 4.2. We have now pi?ip = pm- Since p'^^ is not constant, one hasp|^y < Pm- 

(2) (3) 

Arguments as in Step 4.1 show that we must have p^Jj: — p^Jj: and they are 
attained at the same point Q 7^ O on the contact line, and a contradiction similar 
to (|5.16l) then follows. 

This finishes the proof of Theorem 12.61 
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